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Abstract. Path-integral methods can be used to derive a 'path-decomposition 
expansion' for the temperature Green function of a magnetized free-electron gas 
confined by a hard wall. With the help of this expansion the asymptotic behaviour 
of the profiles for the excess particle density and the electric current density far 
from the edge is determined for arbitrary values of the magnetic field strength. 
The asymptotics are found to depend sensitively on the degree of degeneracy. 
For a non-degenerate electron gas the asymptotic profiles are essentially Gaussian 
(albeit modulated by a Bessel function), on a length scale that is a function of 
the magnetic field strength and the temperature. For a completely degenerate 
electron gas the asymptotic behaviour is again proportional to a Gaussian, with a 
scale that is the magnetic length in this case. The prefactors are polynomial and 
logarithmic functions of the distance from the wall, that depend on the number 
of filled Landau levels n. As a consequence, the Gaussian asymptotic decay sets 
in at distances that are large compared to the magnetic length multiplied by yfn. 



PACS numbers: 05.30.Fk,75.20.-9 



Path-decomposition expansion and edge effects 



2 



1. Introduction 

In a magnetized charged-particle system edge effects are of paramount importance. 
An illustration of this fact is furnished by the phenomenon of Landau diamagnetism 
[jjj, which is due to electric currents flowing near the boundaries of the sample. As a 
further example one may mention quantum Hall systems for which the relevance of 
edge effects has amply been shown. 

The analysis of the influence of the boundary on the properties of a quantum 
many-body system is a difficult mathematical problem. Even if the bulk properties of 
the unconfined system are understood, the presence of the edge leads to a boundary- 
value problem that is often difficult to solve analytically. Leaving out the interparticle 
interaction simplifies this problem quite a lot, although even in that case the analysis 
remains complicated. A physical system that is particularly relevant in this context 
is the non-interacting electron gas in a uniform magnetic field in the presence of a 
confining hard wall. In fact, this is the system in which Landau diamagnetism, with 
currents flowing near the edge, can be studied in its purest form. 

Several methods have been devised to analyse edge effects in the confined 
magnetized free-electron gas. At zero temperature one may try to solve the eigenvalue 
problem in terms of distorted Landau levels and determine the edge currents by 
summing the contributions of the lowest-lying eigenfunctions. Even for a simple flat 
geometry this leads to a rather involved mathematical analysis in terms of parabolic 
cylinder functions Q , || . Recently, we studied the profiles of the particle density and 
the electric current density along these lines Q . 

An alternative approach starts by focusing on the high-temperature regime, where 
Maxwell-Boltzmann statistics applies. In that case a convenient tool is furnished by 
the one-particle temperature-dependent Green function. As shown by Balian and 
Bloch H the Green function for the confined system can be related to that of the 
corresponding system without boundaries by making a systematic expansion that 
accounts for an increasing number of reflections of the particles against the confining 
wall. The ensuing multiple-reflection expansion was used in recent years @, to 
determine the profiles of the particle density and the (electric) current density for 
small values of the magnetic field. These small-field profiles had been found before 
from perturbation theory || , || . It turns out to be difficult to generalize these results 
for the profiles to arbitrary field strength and to relate them to those obtained by 
means of the eigenvalue method. 

Some time ago Auerbach and Kivelson jioj invented a path-integral method to 
analyse boundary effects in Green functions. By suitably decomposing the relevant 
paths near the edge they derived a so-called 'path-decomposition expansion' (PDX) 
for the one-particle Green function. The aim of the present paper is to see whether 
the use of PDX may shed light on the difficulties mentioned above and whether it may 
lead to new results on the profiles of physical quantities for arbitrary field strength, 
both for the high-temperature region and in the regime of high degeneracy. 

The plan of the paper is as follows. We start by a review of the path-decomposition 
expansion and its derivation from the Feynman-Kac path integral. Particular attention 
will be given to the convergence of the PDX series. It will be shown that a suitable 
resummation can greatly enhance that convergence. The connection with the multiple- 
reflection expansion will be established. Subsequently, the extension of the method so 
as to include magnetic fields will be discussed by starting from the Feynman-Kac-Ito 
representation. 
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For the specific case of a non-interacting charged-particle system in a uniform 
magnetic field, and confined by a hard wall parallel to the field, the general form of 
the terms in the PDX series can be established in detail. That result will be used 
to determine the first few terms of the asymptotic expansion for the profiles of the 
particle density and the current density. This asymptotic expansion is valid far from 
the edge and in the high-temperature regime. In contrast to earlier work we will not 
need to restrict ourselves to small field strengths, as we shall establish the full field 
dependence of the profiles. As it turns out, the precise knowledge of the asymptotic 
profiles for high temperatures and arbitrary fields is essential in determining how the 
profiles for the degenerate case depend on the filling of the Landau levels. 



2. Path-decomposition expansion 

Consider a particle in an external potential V(r), i.e. with the Hamiltonian 

H=£ + V(r) (1) 

where we have chosen units in such a way that the particle mass drops out. The 
equilibrium quantum statistical properties of a set of particles moving in the potential 
V is governed by the temperature Green function Gp(r',r), with j3 the inverse 
temperature. Its path-integral representation is given by the Feynman-Kac formula 

G p (r>, r) = (r'\e-P H \r) = J d^f(co) exp - jf" dr V(lu(t)) (2) 

where u>(t) describes the path and d//£ f is the conditional Wiener measure jnj. 

If a hard wall confines the particles to a region of space, the potential can be 
written as V(r) = Vo(r) + V w (r), where V w is a steep wall potential and Vq is a 
smooth external potential. Exact evaluation of (^|) for the confined problem is in 
general not possible, even if the corresponding unconfined problem can be solved 
completely. In this section we will explore the use of the so-called 'path-decomposition 
expansion' (PDX), first introduced by Auerbach and Kivelson |ltj| , to determine the 
Green function of the confined problem. 

To simplify matters, consider the one-dimensional case, with a hard wall at x = 0, 
i.e. V w (x) = oo for x < and V w (x) =0 for x > 0. As a first step, we split the Green 
function into two parts 

G p (x',x) = G° (x',x) + G (x',x). (3) 

Here Gjg is the Green function for the problem without a wall. In order to calculate 
it, one needs to specify the potential Vo(x) for x < as well. We will take the latter 
to be the analytical continuation of the potential for x > 0. We shall assume that the 
resulting Vq{x) is such that G can be evaluated in closed form. 

The second term of (||) is the difficult part. It is a correction that contains 
contributions from all paths crossing the boundary at least once, with an additional 
minus sign so as to compensate the corresponding contributions in G° . In order to 
calculate Gp, one discretizes the path integral in the usual way by introducing n evenly 
spaced grid points at r m = me„, with e n = (3/{n+ 1). Subsequently, one decomposes 
the paths at the boundary pTC( ] . Here 'decomposing' means that the paths are split 
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into two at the point r, where they cross the boundary for the last time. Choosing 
this 'point of no return' between r m and r m+1 , one writes the path integral for G° as 



n />0 poo 

G c Jx',x) = - lim V* / dx m / dx m+1 

H n^oa * — * / / n 

m =l J -°° J ° 



)G° en (x m+u x m )G° men ( ) (4) 

for x and x' both positive. 

In the small interval between r m and r m +i the potential Vq can be ignored, since 
the error will vanish in the continuum limit. Hence, we may use in that interval 
the 'free' propagator G° p(x',x) = (2tt0)~ 1 ^ 2 exp[— [x 1 — x) 2 /20\, where we have put 
h = 1. The free propagator satisfies the identity 

G%(\x'\, -\x\) = sgn(x') lim jf" dr -^jG° f !0 _ T (x', x") G?, r (0, x). (5) 

This identity, which is a generalization of that used in follows directly by 

differentiation of the relation 

f dr[r(l - T )]- 1 / 2 e- a2 / T - b2 / ( - 1 ^ = n Erfc(a + b) (a > 0, b > 0) (6) 
Jo 

with respect to fe. Using (Q) with x < and a;' > in (0), we get 
G^fir/.a;) = — lim lim > / dr / dx m 

x ~q^jjG p-T-rnt^ {% i x )G T (0 , x m )G mtn (x m , x) . (7) 

The integral over x m can be extended to the interval [— oo,oo], if a compensating 
factor 1/2 is inserted. In fact, only small values of x m contribute anyway, at least 
in the continuum limit, owing to the presence of the second G° function. For these 
small values of x m the integrand is approximately invariant under a change of sign of 
x m . Subsequently, we may join the two G° into one, so that we get a closed integral 
relation 

G c (x', x)=- lim ~ J* dr A G/3 _ r(a /, x ") G ^ x ) (8) 

for positive x and x' . This integral relation is the PDX formula derived in JhJ. Since 
the right-hand side contains the original Green function G, we can iterate this integral 
equation by inserting (|^). In this way we arrive at the PDX series: 

Ge(x\x) = G%(x\x)-lim l [\t^-G%_ t {x',x'')G° t {Q,x) 



10 2 J Q dx h 

fl\ 2 f d f ~ r 
+ lim - / dr — lim / dr' 

x"io\2J J dx"x"'ioJ 

x ^ 7 G°_ T _ T ,(x', x'") G° T ,(0, x") G° T (0, x)-... (9) 

To study the convergence of the PDX series we look at the special case of a 
vanishing external potential Vq. The Green function for a free particle in the presence 
of a hard wall can be calculated using a reflection principle. It reads 

GfAx', x) = G%(x', x) - G%(-x', x) (10) 
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for x and x' positive. Since the first term is Gjg, the second term must be the correction 
Gp. To check the validity of (^|) for the present case, we employ (||) repeatedly, with 
the result 



G fil3 (x',x)=G%(x', X )- V2-" )G%(-x',x). (11) 




This is indeed identical to (jlO|). It is clear that all terms in (g) are necessary to 
reproduce the correct result. In addition, we cannot change the order of integration 
and taking the limit in (||). In fact, since one has \im x "io d x >>G°j T , (0, x") — 0, this 
would give an incorrect result. This suggests that (|^) is not the most convenient form 
to use. 

A slightly different series is obtained by modifying (|§|) as follows: 

G c (x',x)=- lim [ P dr^jjG^rix'^'^G^x) (12) 
x T lo J q ux 

where the limit x || is the average of the limits x \ and x J, 0. Since G T (x',x) 
vanishes for x < 0, at least for a hard wall, we have merely added zero to the right- 
hand side of (||). If we iterate (|l2|), with (||) inserted, we get the resummed PDX 
series 

G p (x>, x) = G%{x', x) - Jim J dr A G ° _ T (x', x") G° T (0, x) 

+ lim / dr — — - lim / dr' 

^"IToJo dx"x'"iioJ 

9 -GO_ T _ T ,(x',x"')G T ,(0,x")G° r (0,a ; )-... (13) 



dx'" 

Let us consider again the case Vo(x) = 0. One easily verifies that the correction 
Gp is given by the second term of ( |l3| ) alone. The convergence of the resummed PDX 
series is thus found to be much better than that of of the original one. All higher-order 
terms in ( |l3| ) vanish separately in the present case, since one may prove 

lim A lim fdr'-^GO _ T ,(x' ,x"') G° f T ,(0,x") = 0. (14) 
x"no dx" x"'no J dx"' J,T 

Note that here we are allowed to interchange the order of integration and taking the 
limit. This property is an additional advantage of the series in (|l|). Retur ning to 
the general case with V (x) ^ 0, we expect that both favourable properties of the 
resummed PDX series (fast convergence and invariance under interchange of the order 
of integration and taking the limit) are conserved. Of course, in general the series will 
no longer terminate after the second term. Nevertheless, in some applications only a 
few terms in the expansion are relevant, as we shall see in the following. 
The resummed PDX series (|l^) is of the general form 

oo 

G p (x',x) = Y / G i ;\x',x) (15) 

where we put G^ = G^. The term of order n involves n positions at the boundary. It 
can be seen as arising from paths along which the particle hits the boundary n times. 
These multiple reflections at the boundary form the basis of the multiple-reflection 
expansion, which was derived by Balian and Bloch p] quite some time before the 
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path-decomposition expansion was written down. A close inspection shows that the 
two expansions are completely equivalent. 

Note that in principle the PDX formula (g) and the PDX series (^) can be applied 
to any problem involving distinct spatial regions, for example to tunneling problems 
Jici| |. In contrast, the modified PDX formula ( |l2| ) depends on the presence of a hard 
wall. The application of the resummed PDX series (|l^) is likewise limited to hard- wall 
problems only. 



3. Magnetic Field 



We will now apply the methods of the previous section to a confined free-electron gas 
in a uniform magnetic field. The Hamiltonian is given by 

1 



H=-{p-AY + V w {r) 



(16) 



where V w is again a hard- wall potential. Because of the symmetry of the problem we 
will choose the Landau gauge A = (0, Bx, 0). The factor e/c, with e the charge of the 
particles, has been absorbed in the constant B. 

The presence of the vector potential complicates matters. The Feynman-Kac 
representation (||) of the path integral is no longer valid. We have to use the Feynman- 
Kac-Ito formula instead, which in the special case of V • A — reads [nT| 



rP 

AtV w {u{t)) +i / drw(r) • A(w(t)) 
o Jo 



(17) 



If we let u x denote the ^-component of the path, we can replace the exponential factor 
containing V w by #(inf r oj x {t)). Since the factor that contains the vector potential is 
independent of z, the integral over the z-component of the path gives a trivial factor 
(27r/3) -1 / 2 exp[— (z' — z) 2 /20\. The part of the Green function that depends on x and 
y will be denoted by G ±.p(r' ,r) in the following. 

The path integral over the y-component of the path can be evaluated by a 
Fourier-transform technique. In fact, discretizing the x- and y-components of the 
path, with n intermediate points, we write w(r m ) = r rn = (x m ,y m ), with ro = r and 
i" n +i = r'. The integral in the exponent of ( |l7j ) is then given by dr lu(t)- A(u(t)) = 
Em=o( r ' m + 1 



n+1 



G± ) p(r',r)= Yl / dr m (27re„) _1 exp 



(in Ito's convention). We get 

(^*m Tm— l) 



2c„ 



9(x r> 



x exp [i(y m - y m -i)Bx m -i] S(r n+ i - r'). (18) 

The integrals over y m can now be carried out by using the standard Fourier 
representation of the Dirac 5-function. This introduces an integral over an additional 
variable k. Going back to the continuum limit for the path integral over the x- 
component of the path we arrive at 



with 

Gp{x', x, k) = 



G±,p(/,r) = (27T)- 1 / dfc e ik ( y '- yS >G (x' -k/B, x -k/B,k) (19) 

J — OO 

dn x x 'f(u; x )0{mfuj x (T) + k/ B) C -^ dT ^*^ 2 . (20) 
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This function Gp is the propagator for a particle in a one-dimensional harmonic 
potential with a wall at the position —k/B. 

We are now in a position to use the PDX techniques from the previous section. 
The leading term in the PDX series is found by omitting the wall. In that case the 
propagator Gp becomes jDj: 



G (x',x,k) = 



D 



-l 1/2 



2tt smh((3B) 



exp 



Bx' 



B{x' 2 + x 2 ) 

2tanh(/3S) sinh(/3B) 



(21) 



As a matter of fact, Gp is independent of k, since the only fc-dependence in (|2p| ) is in 
the position of the wall. After performing the integral over k, which is Gaussian, we 
find that the leading term in the PDX series is given by 



B 



4tt sinh(/35/2) 



exp 



B 



r) 2 + —(x' 



x){y' - y) 



(22) 



4tanh(/3B/2) 

which is indeed the Green function in the Landau gauge for the unconfined system. 

The next term in the resummed PDX series ( |T^ ) (or (|l5])) is more complicated. 
The integral over k is again Gaussian (in fact it is Gaussian for all terms), but the 
additional integral over r is not. If we set t\ = tanh(ri?/2), si = sinh(ri?/2), 
ti = tanh((/3 — t)B/2) and S2 — sinh((/3 — t)B/2), we can write 

B 2 rP 



Gi 1) a (r',r) = - 



with 
and 



'-L./3 



16tt 3 /2 



dT fp T ( r 'i r ) eX P 9p\r( r '' r ) 



a/2 



B 

T 



SlS 2 {tl 

x'h+xtz + iiy' -y)] 2 



i 2 ) 1 / 2 



x 

h 



X 

h 



- y) 



h + 



t 2 + 



(23) 



(24) 



(25) 



^>' r) 4 L 

Similar expressions can be found in |Q. Note that the formulas in |?J slightly differ 
from those given above. We have made use of the property Gp(r',r) = [Gp(r,r')]* 
and of the possibility to change r into (3 — r to write f^ and in a form that is 
more symmetric. 

The higher-order terms in the resummed PDX series can be found along similar 
lines. For the special case r' = r they have been collected in Appendix A. They are 
found to agree with those derived in |?J , after appropriate symmctrization. 



4. Asymptotics (non-degenerate case) 

The particle density and the (electric) current density can both be found from the 
Green function. In the absence of quantum degeneracy the particle density is directly 
related to G±.p via 

Po 



p(x) 



-G±jj(r,r) 



(26) 



where po is the bulk density and Z± = B/[4ir sinh(/3£?/2)] is the transverse one-particle 
partition function per unit area for the bulk. The expression for the current density 
is slightly more complicated, involving derivatives of Gj_ j( g: 



^G x ,p(r',r)--H-G^(r',r) 



-Bxp{x). (27) 
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Using only the n = term of the PDX series in the expression for p(x) yields 
the bulk density p , as it should, since ^(r, r) = Zj_. Therefore we will consider 
the excess particle density Sp(x) = p(x) — p instead of p[x) in the following. Since 
there is no bulk current, the n = term of the PDX series does not contribute to the 
current density. 

To determine the profiles of the excess particle density and the current density for 
arbitrary distances from the wall we need to evaluate all terms in the resummed PDX 
series. However, the r-integral in ( |2"5| ) cannot be carried out analytical ly. Likewise, 
evaluation of the multiple r-integrals in the higher-order terms given in |Appcndix A 
is in general not possible. 

For large distances from the wall (in units of the magnetic length 1/yB) the 
leading contribution to the profiles comes from the n = 1 term in the resummed PDX 
series, as we will discuss presently. Moreover, the r-integral in ( p3| ) can be evaluated 
analytically in that limit. It is thus possible to derive asymptotic expressions for the 
profiles of the excess particle density and the current density that are valid far large 
VBx. 

A change of variables y — z 2 /(zq 2 — z 2 ), with z = tanh[_B(2r — /3)/4] and 
zq = tanh(B/3/4), brings (|23|) into the form 



B 3 l 2 x 1 



8^3/2 z y 2 

(1 



■ exp 



x exp 



zo 2 )y 



' Bx 2 ' 




2z 


L 







Ay 



vw + v) 



zo 2 )y 



2z 



Bx A 



(28) 



Because of the presence of Bx 2 in the exponent, only small values of (1 — z 2 )y/z 



contribute to the integral for large \JBx. Since one has < zq < 1, this implies small 
values of (1 — ZQ 2 )y. Note that this does not necessarily mean that y itself is small, as 
zo may be close to 1. For large \[Bx the factor yj\ + (1 — Zo 2 )y in the integrand can 
be replaced by 1. Subsequently, we can use 

dy 



-ay = e »/2K Q ( a /2) 



(29) 



where Kq is the modified Besscl function of the second kind. In this way we arrive at 
the following asymptotic expression for the transverse part of the Green function for 
large \/~B x: 



G% (r,r)«- 



B 3 / 2 x 1-zo 2 

8V27T 3 / 2 zo 3 / 2 



■ exp 



zo 



4:Zn 



Bx 2 



Ko 



The next term in the asymptotic expansion of G^g is 



l-zo 2 



4zn 



Bx' 



(30) 



B 3 l 2 x {\-zp 2 f 
8V2tt 3 / 2 4zo 3 / 2 
1- z 2 



■ exp 



1 



zo 



4zn 



Bx 2 - Ko 



Az 
1-zq 2 
Az 



,0 



Bx' 



Bx' 



(31) 



(1 



zo 2 )y 



1 + (1 -_£o 2 )y/2. Since 
is indeed of 



which can be derived by substituting yl 

z[Kx(z) — K (z)]/ K (z) is bounded for all positive z, we see that ( |31 
higher order in l/{-/B x). 

Having investigated the n = 1 term in the resummed PDX series, we may turn 
to the higher orders. From a detailed analysis (see Appendix A) it is found that all 



Path-decomposition expansion and edge effects 

0-2 r , 



0.1 



0.0 



Figure 1. Numerical results for G±g/G±\ as a function of £ = \fBx, for 



/3B = 4. 



terms with n > 1 are of higher order in in comparison with (|30|). In figure [l] 

we have plotted G^/G^ as a function of £ = \f~B x, for a representative value of 
(3B. The decay is in good agreement with (A.1C). 

The asymptotic expression for the excess particle density at large \[Bx follows 
by substituting (|30|) in (p6h: 



Sp(x) 



B l / 2 x 

-p0 — 7^= z 



1/2 



exp 



1 + V 



Bx 2 K 



1 



Bx' 



(32) 



where we have used that Zj_ is given by Z± = B(l — zo 2 )/(8ttzo) in terms of zq. 
In a similar way an asymptotic expression for the current density at large VB x can 
be derived. In leading order it is found to be proportional to the asymptotic excess 
particle density: 



3y(.x) 



^Bx 5p{x) 



(33) 



with Sp(x) given in ([32|). This simple proportionality ceases to be valid, if higher- 
order terms are incorporated in the asymptotic expansion. Comparing ( |33| ) to ( p7[ ) 
we see that there is a compensation between the term proportional to p{x) and the 
term that contains the derivatives of the Green function. For the n = contribution 
this compensation is complete, but for n — 1 only half of the second term in (|27]) is 
cancelled, at least in leading order in l/{\fB x) . 

It must be stressed that both (|3|) and (|3^) are valid for large \f~B~ x, whereas (3B 
may take arbitrary values. If, apart from Bx 2 , also [(1 — zq 2 ) / zq]Bx 2 is large, we can 
simplify ( |3C| ) to 

~z^ ( Bx 2 ' 



B VT 



exp 



(34) 



z \ 2z 

by using the asymptotic expansion for the modified Bessel function. In that case the 
excess particle density profile is asymptotically given by 

Bx 2 



Sp(x) ~ —po cosh(/3B/4) exp 



2tanh(/3S/4) 



(35) 




Large [(1 — zq 2 )/ zq]Bx 2 implies that the regime zq — > 1 or (3B — > oo is not included, 
whereas no such limitation is imposed on the use of (32). For fixed B this is not a 
serious limitation in the present context of a non-degenerate electron gas, since for 
(3 — > oo we have to use Fermi-Dirac statistics anyway. In the next section it will be 
shown that the Green function in the form ([30|) is crucial to obtain information on the 
asymptotic profiles for a degenerate electron gas with Fermi-Dirac statistics. 

To investigate the validity of the asymptotic expressions for mentioned 
above we have compared ( ^o|) and (34) with numerical results based on (|2^) or (g8j). 
The results are drawn in figure ||. It is clear that for (3B = 4 both asymptotic 
expressions are adequate, even for relatively small values of £, whereas for (5B = 16 
the performance of ( |3C| ) is much better than that of (|34|). 

Finally, for small B, the expression (pq) for the excess particle density yields 

p(x) « -po (l + ^B 2 - e~ 2 ^. (36) 

Indeed these terms correspond to the leading terms (for large x 2 / 0) in the expression 
for the density up to order B 2 , as given in 0. 

In closing this section on the non-degenerate electron gas we remark that the path- 
integral representation can be used to derive a strict bound on the particle density for 
all values of x. Upon setting k' = x — k/B we find from ©, © and ©: 

Sp(x) = r d fc' / d4i(^k'-x-iidu; x (r))e-i B2 ^M^(r)}\ (37) 



All paths that contribute to this integral must pass below the point k' — x, while 
starting and finishing at k' . Now look at paths that go via a point below k' — x precisely 
at t = (3/2. Since these form a subclass of all allowed paths, the corresponding path 
integral provides a lower bound on |5p(a;)|: 

i ) />OC pk —x p 
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iB 2 /| /2 dr[^(r)] 2 



(38) 



The path integrals in this expression are now unrestricted, so that they are given by 
Gp (see (pi])). Integration over k' and x' — k' (in that order) gives 

1 1/2\ 



|5/o(a?) | > ipo Erfc ( a; 



2tanh/3S/4 



(39) 



which is the bound for all x 
implies 



that we set out to derive. In the limit of large x this 



lim 



x exp 



Bx 2 



2tanh(/3B/4) 



\5 P {x)\ > 



Po 



2tanh(/35/4) 



B 



1/2 



(40) 



This inequality is consistent with (|35|), as it should be. In particular, the Gaussian 
decay of Sp(x), with the same characteristic length as in (|35|), is corroborated. 

As can be seen from the results (32) and (33) the decay towards the bulk value 
of both the excess particle density and the current density is Gaussian, modulated 
by a Bessel function and an algebraic factor. For not too large (3B the decay of the 
excess particle density, as given by (^), is strictly Gaussian far from the edge. The 
asymptotic decay of the current density is likewise Gaussian, albeit with an extra 
algebraic factor. The characteristic length on which the Gaussian decay manifests 
itself is proportional to [tanh(/3B/4)/i?] 1 / 2 . As we have shown above, the Gaussian 
decay for the excess particle density is consistent with a lower bound that can be 
derived exactly. For the current density it is consistent with the upper bound on 
the absolute value of the current density that has been derived by Macris et al. |Q . 
However, it should be remarked that the upper bound obtained in that paper is 
rather wide. In fact, the characteristic length of the Gaussian function in their upper 
bound is the thermal wave length, which is independent of the magnetic field. This 
characteristic length is larger than that in the Gaussian found here, at least for non- 
vanishing magnetic fields. 



5. Asymptotics (degenerate case) 

The particle density p{x) FD of a degenerate Fermi-Dirac system at temperature T = 
and chemical potential p is related to the density of the non-degenerate system by a 
Laplace transformation: 



(41) 



' 2Z 
d/ie"*p(a;) fD = —zp(x). 
io PoP 

Here Z = Z\\Z± is the total one-particle partition function per unit volume for the 

bulk, with Zn = (27T/3) -1 / 2 ; the factor 2 takes the spin degeneracy into account. The 



relation (|l|) implies that we can calculate the excess particle density Sp(a 
Sp(x) via an inverse Laplace transform p3| : 

W*)fd = 7T- / d/3 e^— Sp(x) 



'FD 



PoP 



fr( 



(42) 



with arbitrary c > 0. Hence, the asymptotic behaviour of the excess particle density 
of the degenerate system for large \f~Bx can be obtained on the basis of the results of 
the previous section. 
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Let £ = v B x and v = fi/B, and introduce a new integration variable t by writing 
f3 = £(it + !)/-£>, so that c is given by £/-B. If we express the right-hand side of ([32]) 
in the variables £ and and substitute it into we get 



£p(x) 



£3/2^1/2 
167T 3 



df 



i- V 



x exp 



4z 



K 



4z Q 



(it + l) 3 / 2 Z 3 / 2 



(43) 



In the new variables we have z$ = tanh[£;(it + l)/4], so that large £ implies zq ss 1 and 
1 — zq 2 ~ 4exp[— £(it + l)/2]. Consequently, the argument of iio in ( ]l3| ) is small in 



absolute value, so that we can use the series representation 



n=0 



™ 1 Z 

E--7-log(| 



m— 1 



i 



2 2 ™(n!) 2 ' 



(44) 



for the modified Bessel function. In this way we get 



5p(x] 



FD 



-e 



47T 3 



/2' 



4n 



(i* + 
Upon using the identity Jl4| 

e (it+l)x 



n 

E 

m— 1 



l 

m 



dt- 



7 - log 



2-kx u ~ 



n=0 Z 



dt 



(nl) 2 

e [„-(n+l/2)]£(it+l) 
(it + l) 3 / 2 

(i/>0) 



(it + rr v y r(i/) 

we arrive at the asymptotic expression for the excess particle density 

1/2 



(45) 



(46) 



B 3 / 2 £ 
-e 



r 5/2 



["-1/2] 
2 /2 V 

2 2 "(n!) 2 













-i)] 









4[!/-(n+i)] 



n , 

+ E 1 



7 - log 



(47) 



with [x] the largest integer less than or equal to x. The asymptotic expression derived 
here is valid for large £ and fixed v. 

The profile of the current density for large £ and fixed v can likewise be obtained 
from the results of the previous section. In fact, because of the linearity of the inverse 
Laplace transform the asymptotic form of the current density is related to that of the 
excess particle density in the same way as in ([33]): 

j y {x) FD K -]-Bx 5p{x) FD . (48) 



The expressions (47) and ( J48| ) for the asymptotic profiles of the excess particle 
density and the current density are identical to the leading terms of the asymptotic 
expansions in [Q, which have been obtained by solving the eigenvalue problem 
and analysing the asymptotics of the eigenfunctions. It is also possible to recover 
the higher-order terms of |4|] by inserting higher-order terms in the approximate 
expressions for the factors zq and 1 — z 2 in the integrand in (^), taking into account 
corrections like (|3l[), and including more terms in the resummed PDX series as well. 
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The asymptotic behaviour of p7| ) (and of (48)) is Gaussian in £, so that the 
characteristic length is the magnetic length 1 / \f~B for a completely degenerate electron 
gas. Furthermore, the Gaussian is multiplied by a prefactor that depends algebraically 
and logarithmically on £. For v just above a half-odd integer, that is, for chemical 
potentials /i slightly above a Landau level, the profile of the excess particle density 
shows a singular behaviour that is a remnant of the de Haas- van Alphen effect. A 
numerical assessment of the convergence of this asymptotic expression can be found 
in§. 

The dominant term in the asymptotic behaviour comes from the highest Landau 
level with the label \v — 1/2]. Since the prefactor of the Gaussian in this term is 
proportional to £ i l v -l/ 2 ]+ 1 i the onset of the Gaussian decay shifts to larger and larger 
values of £, if v increases. In fact, ( |47|) is useful only for £ 2 v, or equivalently for x 
large compared to the cyclotron radius sfjl/B of particles at the Fermi level. If v is 
large, a different behaviour can be expected in the regime £ 2 « v, before the ultimate 
Gaussian decay sets in at £ 2 ^> v. 

In conclusion, we have studied the edge effects in the excess particle density and 
the current density of a magnetized free-electron gas, which is confined by a hard 
wall. In particular, we have studied the long-range influence of the wall on these 
quantities, by determining their asymptotic profiles both for the non-degenerate case 
and for strong degeneracy. New results have been obtained for both these cases. In the 
former case the asymptotic spatial profiles were found to be Gaussian (or Gaussian 
modulated by a Bessel function), with a characteristic length that is proportional 
to [tanh(/3i?/4)/_B] 1 / 2 . In the latter case the asymptotic behaviour depends on the 
number of filled Landau levels n = [/i/B—1/2]. In fact, it is determined by a Gaussian, 
with a characteristic length equal to the magnetic length multiplied by a 

polynomial and a logarithmic prefactor. Since the degree of the polynomial prefactor 
grows with n, the Gaussian character of the asymptotics comes to the fore only for 
distances that are large compared to ^/ri times the magnetic length. 
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Appendix A. Higher-order terms in the PDX series 

In this appendix we study the asymptotic behaviour of the terms with n > 1 in the 
resummed PDX series for the Green function, for large values of \J~Bx. The general 
form of the term of order n in the resummed PDX series is 

G%{r',r) = H \ n+ f^ +2)/2 J" M ■ ..j\r n 6(r n+1 ) 

with r„ +1 = P — J27=i T «- The functions and grW can be found in Here we 
collect them for the case r' = r, which is relevant for the particle density. In that case 
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we can symmetrize the expressions in ti and i 2 - As a result they get the form 

-(n+l)/2 , 




(A.2) 



and 



In) 



Bx* 



Si 



sinh(J5ri/2; 



+■ *a) 
and 



1 



1 

^2 



h+t 2 



En+l . 
i=l • 

Pochhammer's symbol a(a 



(A.3) 



with ti = tanh(Br,/2), 
l)-..(a + n-l). 

For large \f~Bx the dominant contribution to the integral comes from the 
integration region for which the factor multiplying Bx 2 in \g^ | is minimal. This is the 
case for t\ = ti = (3/2 and Ti = (with 3 < i < n + 1). Therefore, we introduce on a 
par with zo — tanh(_B/3/4) the new integration variables z+ = tanh[_B(/3 — t\ — T2)/4], 
z_ = tanh[_B(Ti — T2)/4] and for n > 2 also Zj = tanh(i?Ti/2) = i< (i = 3, . . . , n). If the 
integrations are carried out in the order Zj, z + and z_, the allowed intervals of these 
variables are z_ e [— z , zo], 2 + <= [0, (z — |z_|)/(l — z |z_|)], and z, £ [0, 2 z + /(l + z 2 )1, 



with an additional condition on resulting from the 9- function in (A.l 



We now have to rewrite the integrand of (A.l) in terms of z + , z_ and Zj. Let us 
consider small values of z + and z^. The function then gets the form 

rz (i - z_ 2 ) 



(n) 



(r,r) « - 



2 

zo 4 z- 



z z - z_ z 
+ z z z_ 



4z 2 z_ 2 + z 2 + z_ 



-z+ + 



(A.4) 



(zo 2 - z_ 2 ) 2 

From the right-hand side it is seen that it is indeed true that only small values of z+ 
contribute to the integral in (A.l), as Bx 2 is large. In turn this implies that all Zj 
have to be small as well, whereas no condition of smallness is imposed on z_. In 
only the p = terms are relevant for large Bx 2 , since these give the terms with the 
highest power of x. As a consequence we can write / ( -™- > as 

- l/2/i _ 2\ / i „ 2 \ 3 / 2 

(r,r)^(i?x 2 r/ 2 Zo } ' ' 

' 7 ; on— 1/2 



An) 



1/2 
1/2 



2z+ 



i=3 



n 



,1/2 



(A-5) 



again for small values of z + and z;. Finally, the ^-function in (A.l) is equal to 
8{2z + — X^i=3 z ») m the neighbourhood of z + = Zj = 0. 
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Since in the approximation considered here g( n ' does not depend on Zi anymore, 
the integral over these variables can be evaluated easily: 

/ \ / \ 1/2 / \ 



n 

\i=3 



2z+/(l+z + 2 ) 



6 2z , 



_(n-l)/2 



i=3 



i=3 



\i=3 



/2 



r(2*+) 



(3n-5)/2 



2"- 1 r[3(n- l)/2] v "" +/ ^ A ' 6 ^ 

To calculate the integral we have extended the upper limit to oo, since the condition 
z + <C 1 guarantees that only small values of Zj contribute anyway. The subsequent 
integral over z+ gets the following form: 

«o-l*-l 



i-*ol*-l 



dz+ (2z+)( 3 "- 5 )/ 2 

Bx 2 Z 4 Z_ 2 + Z 2 Z_ 4 



4z 2 z_ 2 



z 2 + z_ 2 



Xex P — — / 2 "~ 9T9 ^ z + ■ ( A - 7 ) 

Again we can choose 00 for the upper limit, since only small values of z+ are significant; 
the integral can then be carried out trivially. 

We are left with the integral over z- . Leaving it in its original form we arrive at 

on-7/2 1/2 

G&(r,r) « (-)" _ ' (1 *o 2 ) 



n 3 / 2 (Bx 2 ) n ~ 5 / 2 x 2 



(z 2 -Z_ 2 ) 3 ("- 3 / 2 )(l-Z_ 2 ) 1 /2 



x exp 



' (z 4 z_ 2 + z 2 z_ 4 - 4z 2 z_ 2 + z 2 + z_ 2 ) 3 (»- 1 )/ 2 
Bx 2 z (l-z_ 2 ) 



2 z 2 - z_ 2 



(A.. 



A final transformation of variables, by setting z_ = zoy/y/^/1 + y, leads to the 
following asymptotic expression for Gj^o in the regime of large y/B x: 

G («) fr rl „ (_)n£ _f0 (1 ~ Zp ) 



tt 3 / 2 (Bx 2 )"- 5 / 2 x 2 



v/l + (i-z 2 )y 



" Bx 2 ~ 


j 00 dy 


2z _ 


Jo Vvi 1 + y) 


)y 





[2y 2 (l - z 2 ) 2 + y(l - z 2 )(3 - z 2 ) + l]3(»-i)/a 

(1- z 2 )y 



x exp 



2z 



-Bx z 



(A.9) 



-r(i) 



The expression found here looks very similar to (|28|), which is valid for Gj_ 
As before, we may use the fact that only small values of (1 — z 2 )y contribute to the 



integral for large VB x. As a result, one has the asymptotic relation 

71-1 



___°!Y 

Bx 2 J 



(A.10) 



for large Bx 2 . A similar connection formula holds between the asymptotic forms of 
the excess particles density in various orders: 

3 \ n-i 

(A.ll) 
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again for large Bx 2 . Likewise, one derives for the asymptotic forms of the current 
density in various orders: 



We may draw the conclusion that for large \[B x the n = 1 term in the resummed 
PDX series yields the dominant contribution, both for the excess particle density and 
for the current density. 
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